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Abstract 

Travel-times  of  acoustic  signals  were  measured  between  a  bottom-mounted  source  near 
Oahu  and  four  bottom-mounted  receivers  located  near  Washington,  Oregon,  and  Cali¬ 
fornia  in  1988  and  1989.  This  paper  discusses  the  observed  tidal  signals.  At  three  out 
of  four  receivers,  observed  travel  times  at  M2  and  52  periods  agree  with  predictions 
from  barotropic  tide  models  to  within  ±30°  in  phase  and  a  factor  of  1.6  in  amplitude. 
The  discrepancy  at  the  fourth  receiver  can  be  removed  by  including  predicted  effects  of 
phase-locked  baroclinic  tides  generated  by  seamounts.  "  ■  ■ 

Our  estimates  of  barotropic  M2  tidal  dissipation  by  seamounts  vary  between  2  x  10u> 
and  1  X  1018  erg-s-1.  The  variation  by  two  orders  of  magnitude  is  due  to  uncertainties  in 
the  numbers  and  sizes  of  seamounts.  The  larger  dissipation  (1  X  1018  erg-s-1)  is  the  same 
order  as  previous  estimates  and  amounts  to  4  %  of  the  total  dissipation  at  M2. 
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Chapter  1 


Introduction 


We  measure  travel-times  of  acoustic  signals  between  a  bottom-mounted  source  near  Oahu 
and  four  bottom-mounted  receivers  located  near  Washington,  Oregon,  and  California 
in  1988  and  1989  (Figure  1-1).  The  power-spectra  of  basin-scale  travel-times  exhibit 
prominent  variability  at  tidal  periods  [Spiesberger  et  al.,  1989b],  There  are  two  important 
reasons  for  further  study  of  the  tidal  signals.  First,  numerical  models  for  barotropic 
tides  can  be  validated  by  using  these  new  observations.  Second,  it  is  important  to  detide 
tomography  data  to  better  study  other  ocean  fluctuations  at  similar  periods  [ Spiesberger 
et  al,  1989b], 

The  tides  can  be  separated  into  barotropic  and  baroclinic  modes,  and  their  associated 
influences  on  travel-time  can  be  separated  into  current  and  displacement  effects.  The  tidal 
currents  affect  travel-time  by  alternating  between  alliance  and  opposition  to  the  direction 
of  sound  propagation.  The  rise  and  fall  of  the  water  column  affects  the  travel-time  by 
causing  periodic  vertical  iso-therm  displacements  and  pressure  fluctuations.  Worcester 
et  al.  [1990]  has  shown  that  barotropic  current  effects  can  be  quantified  by  using  recip¬ 
rocal  transmissions;  one-way  transmissions  are  used  in  this  experiment,  so  we  quantify 
barotropic  effects  by  comparing  measured  data  to  barotropic  tide  model  predictions. 

We  find  that  acoustic  travel-times  computed  from  barotropic  tidal  models  are  similar 
to  observations  at  three  out  of  four  receivers.  The  discrepancy  at  the  fourth  receiver 
can  be  understood  by  including  the  effects  of  phase-locked  baroclinic  tides  generated  by 
seamounts  near  the  transmission  path.  Generation  mechanisms  are  important  both  for 
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Figure  1-1:  The  acoustic  tomography  experiments  of  1988  and  1989  consist  of  a  source 
located  near  Kaneohe  Bay,  Oahu,  and  four  bottom  receivers  (R.l  to  R4)  whose  approximate 
positions  are  shown.  The  distances  between  source  and  receivers  are  all  about  4000  km. 
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understanding  the  distribution  of  internal  wave  energy  and  the  energy  budget  of  tidal 
dissipation  [Munk,  1966;  Wunsch,  1975;  Bell,  1975;  Lambeck,  1980;  Hendcrshott ,  1981). 

This  report  closely  follows  Headrick  et  al.  [1990]  and  is  organized  as  follows.  Chapter 
2  discusses  the  experiment  and  the  estimation  of  acoustic  travel-time.  Comparisons  be¬ 
tween  observations  and  predictions  from  barotropic  and  baroclinic  tide  models  are  made 
in  Chapter  3.  Chapter  4  predicts  the  global  dissipation  of  tidal  energy  at  M 2  due  to 
seamounts.  Conclusions  follow  in  Chapter  5. 
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Chapter  2 


Data  Acquisition 


The  travel-time  differences  of  pulse-like  signals  are  measured  at  3  minute  intervals  in 
1988  and  2  minute  intervals  in  1989  between  the  source  at  Oahu  and  the  four  receivers 
(Figure  1-1).  Travel-times  are  intermittently  measured  for  a  total  of  between  four  and 
eight  months  (Figure  2-1).  The  transmission  distances  are  about  4000  km  and  the  transit 
times  are  about  40  minutes. 

Spiesberger  et  al.  [1989a, b]  describes  the  maximum  likelihood  estimator  of  phase 
differences  used  to  estimate  differences  in  travel-time.  The  basic  process  is  outlined  in 
Appendix  D.  With  this  estimator,  the  change  in  acoustic  phase  at  133  Hz,  6(f),  is  related 
to  the  travel-time  difference,  Sr,  by, 


2  n6<f> 
o/cuA t ' 


(2.1) 


where  wc  =  2tt  ■  133  (rad-s_1)  is  the  center  radian  frequency  of  the  transmitted  signal,  u>  is 
the  radian  frequency  of  the  travel-time  oscillation,  and  At  is  the  geophysical  time  interval 
over  which  travel-time  differences  are  measured  (184  s  in  1988,  122  s  in  1989)  [Spiesberger 
et  al.,  1989a]. 
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Figure  2-1:  Shaded  regions  indicate  the  times  when  tomographic  data  were  recorded  at 
each  of  the  receivers. 


12 


Chapter  3 


Data  Analysis 

3.1  Barotropic  Model  for  Acoustic  travel-time 

The  change  in  acoustic  travel-time  due  to  barotropic  tides  is  given  by, 

A/|„(C  I’li)  =?■  —y  f  u(s,t  -  —  )•</.?+  ki  f  ( (s,t  - —)(h  ,  (3.1) 

Cq  UI'r  C0  -'Tr  Co 

where  the  geophysical  time  at  the  receiver  is  t,  the  geodesic  between  the  source  and 
receiver  is  lp,  and  the  mean  sound  speed  is  ca.  The  barotropic  current  and  surface 
displacement  at  location  s  on  the  geodesic  at  time  t  are  denoted  by  u(.s,  t  -  and 
( ( s , t  -  ^ )  respectively,  'the  incremental  on  the  geodesic  is  ds,  and  l{s)  is  the  distance 
from  the  receiver.  The  constant  is  k]  =  3.9  X  10-3  [Spiesbergcr  cl  ai,  1989b],  The  effects 
of  finite  propagation  time  are  taken  into  account  by  the  correction  term  (  — ^).  The  finite 
propagation  delay  amounts  to  a  phase  correction  of  order  ±10°  at  semidiurnal  periods. 

The  first  term  in  (3.1)  contains  the  effect  of  propagating  sound  with  and  against 
the  tidal  current.  The  second  term,  which  amounts  to  a  1  or  2  ms  correction,  contains 
the  combined  effects  of  vertical  iso-therm  displacement  and  pressure  fluctuations  which 
accompany  the  rise  and  fall  of  the  sea  surface.  For  barotropic  tides  the  travel-time  change 
is  about  the  same  for  different  acoustic  multipaths  [Spiesbergcr  ct  al. ,  1989a.]. 

The  elevation,  is  obtained  from  the  Naval  Surface  Weapons  Center  (NSWC)  tide 
model  derived  by  Srhwidcrski  [1978].  Unfortunately,  except  for  A/2,  NSWC  models  of  the 
tidal  current  field  are  unavailable.  We  therefore  use  a  linearized  version  of  Srhwidcrski ’s 


[197X]  Laplace  tidal  equations  (LTE)  to  predict  tidal  currents.  The  equations  are, 


^  =  20  cos  Ov  +  7^ 


Hr  sin  0 


IK)  ~  k,{U, 


(3.2) 


W  =  -2ft  cos  0u  -  [^(/vr/-  /3C)]  -  kdv, 

where  A  and  0  are  longitude  and  colatitude  respectively,  «  and  v  are  east  and  north 
current  speeds  respectively,  r;  is  Newton’s  equilibrium  tide  (which  is  proportional  to  the 
astronomical  tide  generating  potential),  g  is  the  earth’s  gravitational  attraction,  Re  is 
the  radius  of  a  spherical  earth,  0  is  the  Earth’s  angular  rate  of  rotation,  and  kd  is  a 
damping  coefficient.  The  constants  a  and  /3  derived  by  Schwidcrski  [1978]  correct  for  the 
solid-earth  tide,  ocean  tidal  loading  of  the  solid  earth,  and  the  associated  gravity  potential 
perturbations  of  all  tides  and  loadings.  The  solution  to  (3.2)  for  u  and  v  is  outlined  in 
Appendix  A. 

Equation  (3.1)  is  decomposed  into  contributions  from  each  tidal  frequency,  uk,  of 
interest,  and  solutions  are  determined  by  collecting  and  integrating  sinu^-t  and  cosu >kt 
terms  separately.  The  resulting  travel-time  change  model  is, 


A7’bt(<,  Fr) 


K 

Y^ek(^n)cosukt  +  ffc(rR)sin  ukt, 

k=  1 


(3.3) 


where  o^Fr)  and  ^(Tr)  are  the  solutions  to  the  separated  integrals,  and  I(  is  the  total 
number  of  tidal  frequencies  calculated. 


3.2  Barotropic  Model  Comparisons  with  Observed  Tidal 
Signals 

3.2.1  Tidal  Signal  Estimation 

Power-spectra  of  the  differences  in  acoustic-phase  exhibit  prominent  variability  at  tidal 
periods  [Spiesbcrgcr  et  al,  1989a, b]  (Figure  3-1).  For  each  receiver,  tidal  signals  are 
estimated  by  fitting  the  8  most  prominent  tidal  components  to  the  data  (Appendix  A, 

Table  A.l).  Specifically,  we  determine  a*  and  bk;  k  =  1,2,..., 8,  in, 

8 

6<f>(tn, l'n)  =  Y2  a/t(I'R)cosu;fc<n  4-  bfc(I'R)  sin  n=  1,2,...,  TV,  (3.4) 

k=i 
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Figure  3-1:  Power  spectral  density  estimate  for  the  133  Hz  phase-differences  at  receiver 
R2.  This  periodogram  uses  a  Bartlett  data  window  with  42  days  of  data  obtained  in  1988. 
The  tidal  peaks  at  the  A'l,  M2,  and  52  frequencies  have  rms  values  which  correspond  to 
travel-time  change  amplitudes  of  8.6,  23.5,  and  12.3  milliseconds  respectively. 

by  linear  least  squares  [Lawson  and  Hanson ,  1974],  where  the  phase-differences  are  mea¬ 
sured  at  times  <n,  and  the  total  number  of  measurements  equals  N.  Depending  on  the 
receiver,  N  equals  55,000  to  135,000,  so  the  least  squares  problem  is  overdetermined. 
Each  least-squares  fit  to  phase-differences  is  compared  to  barotropic  model  travel-times 
to  determine  how  much  of  the  signal  can  be  attributed  to  barotropic  tidal  effects. 

3.2.2  Comparison  with  Model 

Assuming  all  of  the  (3.4)  signal  is  due  to  barotropic  influences,  each  tidal  component 
is  proportional  to  the  time  derivative  of  the  corresponding  (3.3)  model.  Comparisons 
between  observations  and  the  barotropic  model  are  illustrated  by  the  phasors  of  Figure 
3-2.  Each  phasor  magnitude  corresponds  to  a  data-to-model  signal  amplitude  ratio,  which 
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Frequency 

Rl 

Amp(ms) 

R2 

Amp(ms) 

R3 

Amp(ms) 

Rl 

Amp(ms) 

M2 

26 

22 

15 

7 

52 

1 1 

7 

4 

3 

A'l 

10 

T 

1 

1 

3 

Table  3.1:  Predicted  acoustic  signal  amplitudes  for  barotropic  tides.  The  model  ampli¬ 
tudes,  (ejr  +  f iyl2,  are  determined  from  the  constants  in  (3.3)  and  have  rms  errors  of 
about  2  ms. 


from  (2.1 )  is, 

2^r(a  l  +  bp1/2 


u)CukAt(ci  +  f iy/2' 

and  each  phasor  angle  corresponds  to  a  difference  between  data  and  model  signal  phases, 
which  is  given  by, 


tan  1  -  tan  1 


3.2.3  Analysis 

First,  we  discuss  the  comparison  at  two  semidiurnal  periods  (M2  and  52).  Predicted  tidal 
amplitudes  vary  from  3  to  26  ms  at  our  receivers  (Table  3.1).  At  receivers  Rl,  R2,  and 
R3,  observed  amplitudes  are  0.8  to  1.6  times  predicted  amplitudes  (Figure  3-2).  At  R4 
the  observed  amplitudes  are  2.4  and  3.4  times  predictions.  Predicted  and  observed  phases 
are  within  ±30°  at  Rl,  R2,  and  R3  but  differ  by  up  to  47°  at  R4.  Data  at  Rl,  R2,  and  R3 
agree  with  predictions  from  the  barotropic  model  to  first  order,  but  at  R4,  the  barotropic 
model  is  inadequate.  We  return  to  this  problem  in  the  next  section. 

We  next  discuss  the  comparisons  for  one  semidiurnal  period  (A’l).  The  models  for 
barotropic  currents  at  diurnal  periods  are  inaccurate  within  a  ±10°  band  centered  at  30° 
N  (Appendix  A).  The  transmissions  to  all  receivers  intersect  this  band,  so  comparisons 
between  predictions  and  observations  are  problematic.  Nonetheless,  predictions  at  Rl 
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Figure  3-2:  Comparison  of  phase  and  amplitude  of  measured  and  predicted  tidal  signals 
at  four  receivers  (R1  to  R4)  and  three  frequencies  (M2,  S 2,  A'l).  Predicted  tidal  signals 
are  obtained  from  the  barotropic  model  given  by  (3.1).  Each  vector  angle  corresponds  to 
a  measurement  phase  that  leads  the  predicted  phase  by  the  indicated  amount.  The  unit 
circle  in  the  upper  left  hand  corner  illustrates  the  vector  angle  convention.  Vector  magni¬ 
tude  corresponds  to  the  ratio  of  measured  amplitude  to  predicted  amplitude.  Predicted 
magnitudes  are  listed  in  Table  3.1.  An  exact  match  between  measurement  and  prediction 
is  represented  by  a  vector  of  angle  zero  and  length  one  (tip  just  touching  the  unit  circle  at 
zero  degrees).  Solid  arrows  correspond  to  measurements  of  the  combined  1988  and  1989 
record,  and  the  dashed  arrows  correspond  to  separate  1988  and  1989  data  measurements. 
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and  R‘2  resemble  observations  at  tidal  period  K1  (Figure  3-2).  It  may  be  that  K1  and 
K2  tidal  signals  are  similar  to  predictions  because  their  transmission  paths  travel  smaller 
distances  within  3(F  ±  10“  N  than  paths  to  R3  and  IM  (Figure  1). 


3.3  Baroclinic  Tide  Prediction 


3.3.1  Phase-Stable  Internal  Tides 


Apart  from  generation  regions,  baroclinic  tides  are  not  phase-locked  to  the  ATGF,  and 
their  phases  are  incoherent  at  100  km  distance  (approximately  an  internal  tide  wavelength) 
[Ilendershott, \98\}.  In  generation  regions,  baroclinic  tides  may  be  phase-locked  to  the 
ATGF  [Munk  et  al..  1981].  Open  ocean  evidence  of  this  phase-lock  is  seen  in  a  phase- 
stable  M2  baroclinic  tide  over  a  large  guyot  east  of  Hawaii  [Noble  et  al.,  1988].  The 
transmission  path  to  ltd  passes  within  a  few  hundred  kilometers  of  several  large  seamounts, 
one  of  which  is  indicated  in  Figure  1-1.  We  hypothesize  that  the  semidiurnal  signals  at 
R4  do  not  agree  with  predictions  from  the  barotropic  model  because  of  the  presence  of 
phase-locked  baroclinic  tides  near  the  seamounts. 

The  transmission  paths  to  III,  R2,  and  R3  may  be  near  other  large  seamounts,  but 
because  there  are  no  large  discrepancies  between  the  data  and  barotropic  model  predic¬ 
tions,  it  is  unlikely  that  internal  tides  play  a  zero  order  role  in  determining  tidal  signals 
at  these  receivers. 

Appendix  B  summarizes  a  model  for  the  generation  of  baroclinic  tides  from  the  inter¬ 
action  between  barotropic  currents  and  seamounts.  For  baroclinic  tides,  acoustic  travel- 
times  are  modulated  by  sound-speed  fluctuations  associated  with  vertical  straining  of 
iso-tachs  [Spiesbcrger  et  al.,  1989b].  The  effect  for  the  jth  acoustic  multipath,  Tr^'),  is, 


ATbc(<,FR(j))  =  -5  f  Cim (r,t- 

co  ■'I'rIj) 


l(nMr) 


Dz 


dr , 


(3.5) 


where  the  position  along  the  multipath  is  r,  the  vertical  displacement  of  the  internal  tide 
is  Cinti  and  the  vertical  derivative  of  the  potential  sound-speed  is  . 


18 


M2 

Amplit  ude(  ms) 

Phase(u) 

52 

Ainplitude(ms) 

PhaseO' ) 

Barotropic  Model 

7 

113 

3 

324 

Baroclinic  Model 

20 

117 

5 

265 

Combined  Model 

27 

116 

*7 

286 

Measured  Data 

24  ±  1 

1 14  ±  1 

7  ±  0.1 

279  ±  1 

Table  3.2:  Comparison  of  predicted  and  observed  travel-time  changes  at  A/2  and  S 2  for 
receiver  R  '.  The  uncertainties  on  the  measured  data  are  based  on  the  differences  between 
1988  and  1989  least-square  fits.  Phases  are  based  on  a  sinud  time  dependence  with  t  =  0 
at  0000  Greenwich  mean  time  on  1  January  1989. 

3.3.2  Influence  on  Travel-Time  at  Receiver  R4 

Being  the  southern-most  receiver,  R4  has  a  more  east-west  oriented  ray  path  than  the  other 
receivers.  This  orientation  results  in  a  weaker  influence  from  the  generally  north-south 
oriented  A/2  and  52  barotropic  tidal  currents,  so  the  effects  of  phase-locked  baroclinic 
tides  are  more  apparent.  The  predicted  A/2  baroclinic  tide,  generated  from  one  of  several 
seamounts  near  the  transmission  path  to  R4  (Figure  1-1),  has  an  amplitude  and  wave¬ 
length  of  6.2  m  and  98  km  respectively  (Appendix  B)  and  produces  a  21  ms  amplitude 
travel-time  signal  (Appendix  C). 

The  sum  of  the  predicted  A/2  baroclinic  tidal  signals,  generated  from  all  prominent 
seamounts  along  the  transmission  path,  is  added  to  the  barotropic  prediction.  The 
barotropic  plus  baroclinic  prediction  for  R4  is  almost  identical  to  observations  at  A/2 
and  52  (Figure  3-3  and  Table  3.2).  The  prediction  is  improved  at  A'2  but  not  at  N2 
(Table  3.3).  The  frequency  at  N 2  is  close  to  A/2  (Table  A.l),  and  the  observed  amplitude 
at  N 2  is  small.  We  suspect  that  the  barotropic  plus  baroclinic  prediction  does  not  fit 
observations  at  N2,  because  strong  baroclinic  tides  at  A/2  may  be  smeared  into  N2  by 
dynamic  processes  [ Hendcrshott ,  1981]. 
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Figure  3-3:  The  phasors  on  the  left  are  identical  to  the  combined  1988  and  1989  R4  record 
phasors  shown  in  Figure  3-2.  The  phasors  on  the  right  present  the  same  analysis  using 
the  combined  barotropic  and  baroclinic  models  in  Table  3.2. 


N  2 

A'2 

Amplitude(ms) 

Phase(°) 

Amplitude(ms) 

Phase(°) 

Barotropic  Model 

1 

37 

1 

153 

Ilaroclinic  Model 

6 

315 

3 

117 

Combined  Model 

6 

325 

4 

126 

Measured  Data 

2  ±  1 

331  ±  38 

5  ±  2 

153  ±  18 

Table  3.3:  Same  as  Table  3.2,  except  comparisons  are  at  N 2  and  A'2  frequencies. 
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Chapter  4 


Tidal  Dissipation 


The  effect  of  the  ocean  tidal  bulge  on  the  Earth-Moon  system  is  a  well  known  phenomena 
that  is  discussed  extensively  by  Lambeck  [1980].  The  friction-induced  phase  lag  between 
(  and  T)  produces  a  torque  in  the  Earth-Moon  system  that  is  responsible  for  the  length¬ 
ening  of  the  day  and  recession  of  the  Moon  at  rates  of  about  1.8  ms-century-1  and  3.7 
m-century-1  respectively  [Christodoulidis,  et  al.  1988].  The  applied  torque  transfers  an 
amount  of  orbital  momentum  to  the  Moon  equal  to  the  rotational  momentum  lost  by  the 
Earth  [ Iximbeck ,  1977].  The  net  result  of  this  momentum  transfer  is  an  overall  reduction 
in  system  kinetic  and  potential  energy.  System  energy  is  dissipated  for  tne  most  part 
in  oceanic  tides,  with  some  fraction  (which  has  yet  to  be  quantified)  occurring  through 
conversion  to  baroclinic  tides  [Lambeck,  1980].  In  this  chapter  we  estimate  the  barotropic- 
to-baroclinic  conversion  rate  of  the  M2  tide  at  seamounts. 


4.1  Tidal  Dissipation  at  M2  from  Satellite  Observations 


The  rate  of  energy  dissipation  is, 


<1E 

dt 


(5.86^p-  -  9.58^)10‘,uerg  •  s 
dt  dt 


(4.1) 


where  is  the  total  tidal  acceleration  of  the  Earth  and  ^  is  the  acceleration  of  the 
Moon’s  orbit  in  rad-s-2  [ Lambeck ,  1977].  Recent  estimates  for  M2  accelerations  have 
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been  obtained  from  the  orbital  perturbations  of  17  satellites,  and  have  values  given  by, 


=  -1,15  x  l()-'22rad  -s-'2. 


%  =  —9.73  x  10~24rad  -s”2, 


(•1.2) 


[Christodoulidia,  i  t  ni,  1988],  Substituting  (4.2)  into  (4.1)  yields  2.5  x  101<J  erg-s  1  at  M2 
(about  75  %  of  all  tidal  dissipation). 


4.2  Tidal  Dissipation  by  Seamounts 

We  consider  two  models  for  the  distributions  of  seamounts.  In  the  first  model  the  number 
of  seamounts  with  a  basal  radius  r  >  It  (km)  is  given  by, 

N(It)  =  1,1 1  x  !0,;exp[— 0.G3/?].  (1.3) 

liquation  (4.3)  is  an  extrapolation  of  Jordan  ci  al.  ’s  [1983]  results,  for  the  eastern  Pacific, 
to  the  world  ocean.  In  the  second  model,  the  number  of  seamounts  is  estimated  from 
altimeter  measurements  by  Seasat  [ Craig  and  Sandwell,  1988].  Seasat  detects  many  large 
seamounts  which  are  absent  in  the  model  by  Jordan  el  al.  [1983],  but  it  does  not  detect 
small  seamounts,  as  the  altimeter  profile  sensitivity  to  seamounts  falls  off  rapidly  below  a 
basal  radius  of  20  km  [Craig  c.t  al.,  1988].  These  two  seamount  distributions  are  outlined 
in  Table  4.1. 

The  rate  of  energy  conversion  from  barotropic  to  baroclinic  tides  at  each  seamount  is 
approximately, 

d  y  _ 

-jj-  =  2  (Ax)CgiEi,  (4.4) 

where  Ax  is  the  horizontal  dimension  of  each  seamount  perpendicular  to  the  direction  of 
barotropic  current,  cgi  is  the  group  speed  of  the  internal  tide  (Equation  (B.22)),  and  E,  is 
the  average  energy  density  of  the  internal  tide  (Equation  (B.21)).  The  dimension,  Ax,  is 
obtained  from  the  seamount’s  basal  radius  according  to  Table  4.2.  If  we  take  the  barotropic 
current  speed  to  equal  2  cm  s-1  and  the  inertial  frequency  to  equal  0.73  x  10-4  s-1  (30° 
N),  then  equations  (4.4)  and  (B.23)  yield  a  total  energy  conversion  rate  of  2  X  1016  and 
1  X  1018  erg-s-1  for  the  exponential  and  Seasat  distributions  respectively.  These  estimates 
are  proportional  to  Ax  and  the  squares  of  both  seamount  height  and  the  barotropic 
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Basal  Radius 
r(km) 

Number  of  Seamounts 

Exponential 

Model 

Seasat 

Model 

1  -  2 

360, 000 

0 

2-4 

290,000 

0 

4-6 

83,000 

0 

6-10 

30,000 

200 

10  -  15 

2500 

1200 

15  -  20 

110 

6000 

20  -  26 

5 

8800 

26  -  32 

0 

6800 

32  -  38 

0 

4800 

38  -  44 

0 

2800 

44  -  50 

0 

1600 

Table  4.1:  Seamount  distribution  for  the  world’s  oceans  as  a  function  of  basal  radius  by 
Jordan  et  al.  [1983]  (exponential  model)  and  Craig  and  Sandwdl  [1988]  (Seasat  model). 
The  Seasat  model  is  obtained  by  converting  Craig  and  Sandwell’s  [1988]  histogram  of 
peak-to-trough  separation,  p,  to  basal  radius  using  r=(0.83)p  in  accordance  with  mea¬ 
sured  seamount  slopes.  The  numbers  so  obtained  are  multiplied  by  four  to  account  for 
undetected  seamounts  [Craig  et  al.,  1988].  Some  of  the  discrepancy  between  models  may 
be  due  to  large  seamounts  in  the  western  Pacific  that  were  not  represented  in  the  Jordan 
et  al.  [1983]  survey. 


Rectangular  Seamount  Dimensions 

Exponential 

Seasat 

Dimension 

Model 

Model 

Ax 

(0.62)r 

(0.62)r 

Az 

(0.21)r 

2500  m 

Table  4.2:  Seamount  dimensions  for  the  exponential  and  Seasat  models  (Table  4.1).  The 
length  and  width  of  each  seamount  is  Ax,  and  the  height  is  A z.  Basal  radius  is  denoted 
by  r.  The  scale  factor  for  Ax,  (0.62),  is  derived  from  a  mean  flatness  of  0.31  reported  for 
85  seamounts  in  the  Pacific  [Smith,  1988].  The  Az  scale  factor  for  the  exponential  model 
is  derived  from  the  mean  height-to-basal-radius  ratio  of  0.21  [Smith,  1988].  The  height  is 
set  to  2500  m  for  the  Seasat  model. 
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current  speed  (Equation  (11.23)).  The  exponential  model  yields  insignificant  dissipation 
because  it  contains  relatively  few  wide  seamounts  (Table  4.1).  The  Seasat  model  yields  a 
dissipation  that  is  about  4  %  of  the  total  at  A/2.  Munk  [1966]  and  Bell  [1975]  obtained 
rates  of  5.0  X  1018  and  2.5  X  1018  erg-s-1  respectively  for  the  conversion  to  baroclinic  tides 
by  bottom  roughness.  These  estimates  are  the  same  order  as  we  obtain  using  the  Seasat 
distribution. 
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Chapter  5 


Conclusions 


We  reach  the  following  conclusions: 

1.  Because  barotropic  signal  predictions  at  Rl,  R2,  and  R3  compare  favorably  with 
observations  and  because  barotropic  wave-lengths  have  basin-scales,  the  principle 
discrepancy  between  barotropic  signal  predictions  and  observations  at  R4  is  probably 
not  due  to  errors  in  Schwiderski’s  [1978]  model.  Furthermore,  because  barotropic 
plus  baroclinic  signal  predictions  compare  favorably  with  observations  at  R4,  we 
believe  the  Schwiderski  [1978]  model  is  adequate  for  barotropic  predictions. 

2.  Seamounts  may  be  responsible  for  a  small  (~  5  %),  but  non-negligible,  fraction  of 
the  dissipation  of  barotropic  tides.  The  dissipation  is  sensitive  to  the  heights,  radii, 
and  numbers  of  seamounts  (all  poorly  known). 

3.  The  generation  of  baroclinic  tides  at  seamounts  can  be  accounted  for  with  a  two- 
layer  model  of  the  ocean  (Appendix  B). 

4.  Barotropic  tidal  currents  can  be  predicted  from  Schwiderski’s  [1978]  tidal  elevation 
predictions  with  our  model  (Appendix  A),  except  within  ±10°  of  30°  North  and 
South  (diurnal  frequencies)  and  near  90°  N  (semidiurnal  frequencies). 

5.  Comparison  of  measured  tidal  signals  (from  basin-scale  acoustic  transmissions)  with 
predicted  tidal  signals  (from  barotropic  models)  can  be  used  to  detect  baroclinic 
tides  which  arc  phase-locked  to  the  ATGF. 


(i.  It  is  difficult  to  make  a-priori  prodictions  of  tidal  signals  in  basin-scale  transmissions, 
because  generation  regions  of  baroclinic  tides  are  poorly  mapped. 
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Appendix  A 


Obtaining  Tidal  Currents  from 
Elevations 


Wo  sol vo  (.'{.2)  for  u  and  v  by  defining. 


u  =  b\  cos  cJt  +  b2  sin  u )t, 
v  -  a  1  cos ut  +  c^sinmt, 


( A.l) 


and, 

=9iCos^  +  !,2sin^,  (A2) 

~n  [  -^Oj  =  cos  ut  +  22sinuT, 

where  u>  is  the  tidal  frequency  of  interest  and  gj  ,  02,  b\ ,  62,  y\ ,  j/2 •>  21 >  ar>d  z2  are  constants, 
'the  elevation,  C,  and  Newton’s  equilibrium  tide,  77,  are  periodic  with  frequencies  shown 
in  Table  A.l.  The  y's  and  z's  are  computed  from  (A.2)  by  expanding  known  functions,  C 
and  r/,  in  terms  of  sines  and  cosines.  Substituting  (A.l)  and  (A.2)  in  (3.2)  leads  to  four 
linear  equations,  in  the  unknowns  a\,  a2,  6j,  and  b2,  which  can  be  written  in  matrix  form 
as, 


0  -20  cos  6  -u 

-u ’  kj  0 

—20  cos  0  0  k,i 

kj  u  20  cos 


kd  ^ 

U) 

20  cos  0 

a2 

22 

LJ  ! 

by  : 

y\ 

0  ) 

V  b 2  ) 

V  ‘'1  J 

(A. 3) 
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Semidiurnal  Modes 

Diurnal  Modes 

Component 

Frequency 

(cpd) 

Component 

Fre(}uency 

(cpd) 

M2 

1.93227 

K1 

1 .00273 

S2 

2.00000 

01 

0.92954 

N2 

1.89600 

PI 

0.99726 

I\2 

2.00547 

Q1 

0.89325 

Table  A.l:  Major  Tidal  Frequencies 


The  solution  to  this  matrix  equation  determines  u  and  v. 

We  use  the  values, 

o  =  .69,  13  =  .90,  (A.4) 

derived  by  Schwiderski  [1978],  and  the  value  of  the  linear  damping  coefficient,  kd,  is 
selected  by  minimizing  the  phase  difference  between  predictions  of  AT{t,  M2,Tra)  using 
Schwiderski  [1983]  and  damped  LTE  currents.  The  resulting  best  value  is, 


kd  =  9.0  x  10'V1. 


(A.5) 


The  value  of  k d  also  controls  the  Q  of  the  LTE  model.  We  estimate  Q  by  comparing 
the  equations  of  motion  to  a  forced  linear- damped  oscillator, 


d2x  dx  2 
—  +  k—  +  u^x  =  A  cos ut. 


whose  Q  is  defined  as. 


Q  = 


wr 
k  ’ 


where  the  resonant  frequency  is, 


-rt  =  K2~y)1/2, 


(A. 6) 


(A. 7) 


(A. 8) 


[ Marion ,  1970].  By  letting  east  be  the  r-direction.  and  substituting  (A. 2)  into  (3.2),  we 
obtain  the  following  x-direction  momentum  equation, 


x  Ox 

-7^2  +kd—  -  f  v  =  y}  cose at  +  y2  sin  ut. 


(A. 9) 
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where  /  =  211  ros  0  is  the  inertial  frequency.  Since, 


\%\  %  1.1  x 

\fv\  %  1.-1  x  10“**, 

(y'i  +  3/2) 1/2  ~  4.4  x  10-7, 

\kdu\  %  9.0  x  10~8. 


the  dominant  terms  in  (A. 9)  can  be  written 


du 

dt 


fv  =  0. 


Substituting  the  solution  to  (A.  11), 


and  the  relation, 


ILJU 


u  =  iu>x, 


(A. 10) 


(A.l  1) 


(A. 12) 


(A. 13) 


into  (A9),  we  obtain  the  approximation, 


d2x  dx 

-Qp  +  kd-^  +  oj2x  ~  yi  cos ut  +  t/2  sinwi. 


(A. 14) 


Comparing  (A. 6, 7, 8)  with  (A. 14)  we  see, 

(cj2  - 

Q  «  - - 2-L— 


u  1.41  X  10  4s-1 


kd 


kd  9.0  x  10_6s 


^rr  =  16’ 


(A.  15) 


where  1.41  x  10-4s-1  is  the  frequency  of  the  M2  tide.  This  agrees  with  estimates  ranging 
from  5  to  25  provided  by  llendershott  [1981]. 

A  comparison  of  the  damped  LTE  model  with  the  Schuiiderski  [1983]  model  for  M2 
in  the  North  Pacific  ocean  demonstrates  good  agreement  (Figure  A-l).  Both  models  also 
generally  agree  with  the  data  presented  by  Luyten  and  Stommel  [1990]  extracted  from 
deep  moored  current  meter  records.  By  extension,  we  believe  that  the  damped  LTE 
model  should  be  an  accurate  first  order  predictor  of  the  52,  N 2,  and  K2  semidiurnal 
currents. 

The  damped  LTE  model  does  not  perform  as  well  as  the  Schwiderski  [1983]  M2  current 
model  in  the  vicinity  of  amphidromic  points.  The  finite  difference  approximations  we  use 
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Figure  A- 1 :  The  M2  barotropic  tidal  current  in  the  northeast  Pacific  Ocean  at  0600 
Greenwich  mean  time  on  1  January  1989  as  determined  from  (a)  the  Schwiderski  [1983] 
model  (solid  arrows)  (b)  linear-damped  Laplace  tidal  equations  (Equation  (3.2);  dashed 
arrows)  (c)  data  extracted  from  deep  moored  current  meter  records  (two  arrows  desig¬ 
nated  by  circles  around  their  tails).  The  current  meter  analysis  is  taken  from  Luyten  and 
Stommcl  [1990] 
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for  the  spatial  derivatives  in  (3.2)  are  not  accurate  near  these  areas  of  rapidly  changing 
phase.  These  inaccuracies  are  limited  to  a  100  km  (the  grid  size  of  Schwidorski's  model) 
radius,  and  the  net  effect  on  predicted  travel-time  change  is  less  than  5  %  in  amplitude 
and  3  degrees  in  phase  for  a  -1000  km  geodesic  passing  through  an  amphidrome. 

The  damped  LTE  model  has  an  additional  difficulty  for  diurnal  currents  within  a  ±10° 
band  centered  at  the  inertial  latitude, 


0j  =  sin  *(wfc/2fi). 


(A.  10) 


For  diurnal  tides,  w*.  ~  1  cpd  and  6j  =  30°  N,  and  for  semidiurnal  tides,  u>k  ~  2  cpd  and 
Of  =  90°  N.  The  difficulty  can  be  seen  by  setting  kd  =  0  and  solving  (3.2)  for  u  and  v. 
The  resulting  equations, 


u  = 


g^ncosO^an-PO+^&^-Pa 
li(4Q2  cos2  0—uj2 ) 


PC) 


(A. 17) 


have  denominators  equal  to  zero  at  the  inertial  latitude.  These  zeros  are  canceled  by 
numerator  zeros  when  u,  v,  and  (  are  determined  exactly  from  the  same  equations  of  mo¬ 
tion,  but  a  singularity  occurs  when  the  numerators  are  approximated  by  finite  differences. 
The  finite  fcj  we  use  in  (3.2)  removes  the  singularity  in  (A. 17),  but  the  solution  is  still 
inaccurate  near  30°  N  (Figure  A-2). 
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Figure  A-2:  The  AT  barotropic  tidal  current  in  the  northeast  Pacific  Ocean  at  0600 
Greenwich  mean  time  on  1  January  1989  as  determined  from  the  linear-damped  Laplace 
tidal  equations  (Equation  (A.2)).  Note  the  irregular  performance  of  the  model  in  the 
vicinity  of  30.1°  N. 
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Appendix  B 


Guyot- Generated  Internal  Tides 


In  order  to  determine  the  coupling  between  barotropic  tidal  currents  and  internal  tides 
generated  by  guvots,  we  use  the  two-dimensional  two-laver  model  illustrated  by  Figure  B- 
1.  This  model  is  applicable  to  sea  plateaus  whose  major  horizontal  extent  is  perpendicular 
to  the  direction  of  the  barotropic  current.  We  follow  the  procedure  used  in  deriving 
the  coupling  between  surface  tides  and  coastal  generated  internal  tides  [ Rattray ,  I960]. 
Variables  are  a  function  only  of  the  distance  from  the  center  of  the  plateau  (y-direction) 
in  this  two-dimensional  model. 

A  rigid-lid  approximation  is  used.  The  momentum  equations  are, 


dtv'  +  fu'  =  0, 


(B.l) 


Ot  v"  +  fu"  =  -grdy  C", 


(B.2) 


and. 


dtu'  -  / v'  =  0, 
dtu"  -  fv"  =  0, 


(B.3) 


where  u  is  the  current  parallel  to  the  plateau,  (')  terms  correspond  to  upper  layer  and  (") 
to  lower  layer  terms,  /  is  the  inertial  frequency,  and, 


«7r  =  fifAp/p, 


(B.4) 


is  reduced  gravity,  where, 


Figure  B-l:  Definition  Figure  for  a  two-layer  ocean  model  used  to  predict  the  coupling 
between  a  barotropic  tidal  current  and  the  internal  tide  generated  bv  a  seamount  whose 
extent  perpendicular  to  the  direction  of  the  current  is  considered  infinite. 


Subtracting  (B.2)  from  (B.l),  we  find, 

Ot(v'  -  v")  +  /(«'-  u")  =  grdyC\  (B.6) 


and  the  continuity  equations  for  each  layer  are, 


dy(h'v')  =  QtC\ 

dy{h"v")  =  -dtC- 


(B.7) 


Let  v,  -  (a/d)expi(u>t  +  0),  which  is  essentially  constant  for  any  distance  y  small  com¬ 
pared  to  the  wavelength  of  the  surface  tide.  We  have  contributions  to  the  currents  in  the 
two  layers  from  both  the  barotropic  current  and  the  currents  associated  with  the  internal 
tide.  Thus, 


v'  =  v,  +  v\  u'  =  u,  +  u' 
v"  =  v,  +  v'l  u"  =  us- (- 


(B-8) 


C  =  0, 


c" 

’’lilt  ■ 


(B.9) 


and  from  the  rigid-lid  approximation, 


We  choose  the  following  solutions  which  satisfy  ( 1 J .li )  and  ( H .7 ) , 


y  <  -1 

-f  <  y  <  f 

y  >  f 

>ini 

- B<lk 

A  sin  k j  y 

B(-'hy 

i  M 

1/c/i 

1  /  d2 

v' 

J 

-( u/k2h')Bc'k ™ 

-(iu)/k\ h')A  cos  k\  y 

-( uj/k2h')Be-,k 

r" 

(u>  /  k2h"2)B('k'2V 

(iui/ki  h'[)A  cos  k\  y 

(ui/k2h2 ) Bc~ikiy 

M* 

-if  /u>d2 

-if/ujdi 

-i  f  /u>d2 

«:• 

( if/k2h')Bc,k « 

-(//£,  5'M  cos  fc,2/ 

( if/k2h')Be~ik w 

< 

-( iflk2h"2)Belk 2» 

( // k2h")A  cos  k , y 

-( if/k2h'l)Be-'k «, 

with  the  common  multiplier, 

a  exp'(w,+<^  .  (IU1) 

To  satisfy  ( B.O), 

J  =  ^2-  (B-12) 

In  order  to  patch  the  solutions  at  the  edges  of  the  plateau,  Q"  must  be  continuous  at  y  =  (. 
and  y  =  -C.  Thus, 

Be"*2*  =  Asia  kit.  (B.13) 


The  solutions  (B.10)  already  satisfy  the  conservation  of  total  mass  flux  h'v'  -f  h"v"  - 
dv ,  +  h'v't  +  h"v"  across  y  =  t  and  y  =  —l,  but  conservation  of  p'  mass  flux  requires  that 
h'v'  =  h'( vs  +  v't)  be  continuous  also.  To  satisfy  this  requirement  we  let, 


A  =  A0e~’ia-^2\ 

B  =  B0e -Hv-v/t-kit)' 


(B.14) 


Substituting  (B.14)  into  (B.13)  leads  to, 


Aq  sin  k\i  =  Bo. 

Using  these  relations  and  conservation  of  p'  mass  flux  we  find, 
h'  uj  .  h'  ioj 

- - —  yto( sin  a  +  i  cos  a)  sin  k\l  =  - - —  To(sin  a  +  i  cos  o)  cos  k\t. 

n  2  k2  d  i  k\ 

Bv  equating  the  imaginary  part  of  (B.16)  we  obtain, 

it, 

tan  a  =  ~  tan  fc,t, 
k2 


(13.15) 


(B.16) 


(B.17) 
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\y\<t 

Isrl  >  ( 

Amplitude 

l‘hase(°) 

Amplitude 

Phase(°) 

c" 

lit 

«/lo  sin  k\y 

—  a  +  90 

aB0{\y\/y) 

4>  -  a  +  90  -  k2(\y\  -  l) 

a/di 

a/d2 

4> 

cl<jjAq  cos  k\ylk\h' 

4>  —  a 

au>Bo/k2h' 

<f>  —  a  +  270  -  k2(\y\  —  t) 

wjjAq  cos  k\  yj k\  h'[ 

<t>  —  a  +  180 

auBo/k2h2 

<f>-  a  +  90  -  k2(\y\  -  i) 

a  f  /u>d  i 

d>  +  270 

a  f  ju>d2 

4>  +  270 

a f  Ao  cos  k\y/k\h' 

<fi  —  a  -f  270 

a  fBo/k2h' 

<t>  -  a  +  180  -  k2(\y\  -  £) 

< 

a  f  Aq  cos  k\y/k\h '{ 

</>  —  a  +  90 

afBQ/k2h'2 

<f>  —  ot  —  k2(\y\  -  £) 

Table  13.1:  Generalized  two-layer  internal  wave  solution  for  a  barotropic  north-south  tidal 
current  propagating  over  an  undersea  plateau  of  infinite  east-west  extent. 


and  by  equating  the  real  part,  we  can  write. 


k\  h'  cos  a 
u cos k\t  ’ 


(B.18) 


which  leads  to, 

Z?o  =  (l-l)M!sina.  (B.19) 

«2  « l  OJ 

The  result  of  this  solution  is  a  standing  interna]  wave  over  the  extent  of  the  plateau  which 
patches  to  traveling  internal  waves  at  either  end.  The  generalized  solution  is  collected  in 
Table  13.1,  and  the  particular  solution  for  a  guyot  located  at  30.3°N,  222.5°E  is  summarized 
in  Table  13.2. 

The  cycle-average  baroclinic  energy  density  for  this  two-layer  ocean  is, 


B,  =  f  £  int)V  +  h'U)2  +  h'{v\)2  +  +  h"{v''f)dt,  (B.20) 

where  T  is  the  period  of  the  internal  tide.  Substituting  the  Table  B.l  amplitudes  and 
phases  in  (B.20)  results  in, 


=  ^P9r(v,d2BQ)2 
2  \  u ; 


2  _  p 


(B.21) 


where  the  resulting  energy  density  is  in  joules  per  square  meter.  This  baroclinic  energy 
propagates  at  a  group  speed  that  from  (B.12)  is  given  by, 

u;2-/2 


Cgi  — 


k2u 


(B.22) 
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Constituicnt 

M2 

52 

Amplitude 

Greenwich 

Phase(°) 

Amplitude 

Greenwich 

Phase(°) 

A" 

^int 

6.21  m 

166  +  k2y 

2.08  m 

202  + 

Vs 

2.3  cm-s-1 

68 

0.7  cm-s-1 

104 

1.9  cm-s-1 

346  +  k2y 

0.6  cm-s-1 

22  +  k2y 

< 

0.3  cm-s-1 

166  +  k2y 

0.1  cm-s-1 

202  +  k2y 

us 

1.2  cm-s-1 

158 

0.4  cm-s-1 

194 

< 

1.0  cm-s-1 

76  -)-  k2y 

0.3  cm-s-1 

112  +  k2y 

u“ 

0.2  cm-s-1 

256  +  k2y 

0.1  cm-s-1 

292  -f-  k^y 

Table  13.2:  Two-layer  A/2  and  S 2  internal  tide  solutions  for  regions  north  of  a  22  by  111 
km  undersea  plateau.  Seamount  is  located  at  30.3°N,  222. 5°E  (Figure  1-1),  /  =  11  km, 
d\  —  2200  m,  d2  =  4700  m,  h'  —  700  m,  gr  =  .0059  m-s-2,  A:2(M2)  =  6.41xl0-5  m-1 
(wavelength  is  98  km),  and  k2(S2)  =  6.71xl0-5  m-1  (wavelength  is  94  km). 


The  group  speed,  energy  density,  and  horizontal  extent  of  each  seamount  are  combined 
in  (4.4)  to  obtain  the  barotropic-to-baroclinic  tidal  energy  conversion  rate.  Expansion  of 
(4.4)  in  terms  of  fundamental  parameters  results  in, 

2  /^2  _  /2\ ‘/2 


where, 


f  =  ^,o*w'v(.-f) 

rdxh”\^2 


tan  q 


d2h'{) 


tan 


f2 


/i" 


1/2 


sin2  a. 


(B.23) 


(B.24) 


gTh'h'[  ) 

The  internal  tide  in  Table  B.2  has  an  energy  density  and  group  speed  of  about  1.5  X  105 
erg  •  cm-2  and  160  cm-s-1  respectively.  The  energy  density  near  this  seamount  is  about 
17  %  of  that  contained  in  the  parent  barotropic  tide,  and  the  baroclinic  conversion  rate 
is  about  5.3  X  10H  erg-s-1  (53  Mw). 

The  reduced  gravity  (B.4)  used  in  the  dispersion  relation  (B.12)  is  subject  to  the 
choices  for  the  representative  densities  p'  and  p"  in  the  upper  and  lower  layers.  Values  of 
gT  of  order  10-2  m-s-2  are  typical  for  two-layer  ocean  models,  but  using  gT  =  10-2  m-s-2 
is  not  optimal  for  minimizing  the  mean  square  error  between  the  combined  model  and  the 
measured  travel  time  signals  at  R4.  The  optimal  value  of  reduced  gravity  chosen  for  this 
model,  as  illustrated  by  Figure  B2,  is  about  5.9  x  10-3  m-s-2. 
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Appendix  C 


Baroclinic  Tide  Travel-Time 
Model 


Equation  (3.5)  predicts  the  travel-time  modulation  of  a  single  multipath  caused  by  one 
seamount.  The  multipaths  oscillate  about  the  sound  axis,  and  a  typical  wavelength  (upper 
plus  lower  loop)  is  50  km.  The  potential  sound  speed  gradient  in  the  region  of  the  Table 
B.2  plateau  is  approximately, 


0.13s-1 

0.0s-1 


above  axis 
below  axis, 


(C.l) 


where  positive  z  is  upward.  We  assume  the  average  ray  path  is  above  axis  for  50  %  of  its 
travel  adjacent  to  the  111  km  length  of  the  plateau,  so  the  average  change  in  travel-time 
at  the  M2  period  is, 


ATbc(<) 


_  (6-2 c,uj{t  <o)+0<'n'  m)(0.13s  ])(lllkm)  _  *-(*-£)+««„ 

(2)(  1479  m  s-1)2 


ms, 


(C.2) 


where  #<;int  is  the  phase  of  the  M2  internal  tide,  and  u  is  the  M2  frequency. 

An  evaluation  of  N,  seamount-generated  tides  along  the  geodesic  between  source  and 
receiver  will  produce  N,  baroclinic  signals  at  each  tidal  frequency.  The  vector  sum  of 
baroclinic  models  plus  the  barotropic  model  results  in  the  combined  model, 

N, 

A7’tot(/,^)  =  ATu(t,uk)  +  &Thc(n,t,ujk),  (C.3) 

n  =  l 
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whore  A7|,,-( » ,  t , u.^. )  is  the  baroclinic  signal  prediction  for  tlie  7tth  seamount  at  frequency 
^k- 
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Appendix  D 


Signal  Processing 

D.l  The  Transmitted  Signal 

Pulse  compression  techniques  described  by  Spiesbcrger  et  al.  [1980]  and  Metzger  [1983]  are 
utilized  to  increase  signal-to-noise  ratio  without  degrading  the  phase  difference  resolution. 
The  transmitted  signal  has  a  center  frequency  of  133  Hz  and  a  bandwidth  of  about  17  Hz. 
The  center  frequency  is  phase  modulated  every  8  cycles  using  a  511  digit  maximal  shift 
register  sequence  having  a  period  of  511(8/133)  ss  30.7  s  [Spiesbcrger  el  al.,  1989a].  The 
ideal  transmitted  waveform  is. 

s(t)  so  cos[2x/c/.  -  ^(r*(t)  -  ^)]  (D.l) 

[Spiesbcrger  et  al.,  1980],  where  fc  is  the  center  frequency,  rj;  represents  the  shift  register 
sequence  of  ones  and  zeros,  and, 

22  qi 

*(0  =  £>(*-—),  (D.2) 

where  u (t)  is  the  unit  step  function. 
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D.2  Post-Detection  Processing 

D.2.1  Demodulation  and  Averaging 

The  received  signal,  z(t),  can  be  approximated  by, 

Ntnp 

z(t)  =  ^  bns(t  -  rn),  ( 1 ) - 3 ) 

n=l 

which  is  the  weighted  sum  of  Nmp  multipaths  each  having  gain  bn  and  delay  rn  [Spiesberger 
ct  ai,  1980].  Demodulation  of  the  received  signal  is  achieved  by  sampling  at  4  •  133Hz  = 
532  Hz.  Samples  of  an  ideal  single-path  arrival  are, 

I17T,  7T  117T  7 T 

*[n]  =  cos(--)cos[-(rA..[n]  -  1/2)]  +  sm(  — ) sin[-(rfc[n]  -  1/ -/j,  (D.4) 

where, 

OO 

Ar[n]  =  51  >»[n  -  32/],  (D.5) 

1  =  0 

and  u[n]  is  the  unit  step  sequence.  Because  the  sampling  frequency  is  four  times  the 
carrier  frequency,  the  even  numbered  samples  are  the  cosine  of  the  demodulated  phase 
angle  (weighted  by  ±1)  and  the  odd  numbered  samples  are  the  sine  of  the  demodulated 
phase  angle  (weighted  by  ±1).  The  demodulated  phase  angles  are  low-pass  filtered  by 
block  averaging  to  achieve  a  sequence  of  complex  demodulates  defined  as, 

c{j]  =  x[j]  +  iy[j ];  j  =  0,1,..., J  -  1,  (D.6) 

where  the  square  root  of  minus  one  equals  i, 

3 

x[j)  =  52(-l)m+M2(m  +  4j)],  (D.7) 

m=0 

and, 

3 

y[j\  =  +  4i)  +  J],  (D.8) 

m= 0 

[Spiesberger  ct  nl. ,  1980].  The  resulting  sample  rate  is  four  complex  demodulates  per 
digit,  so  there  are  ./  =  4-511  =  2044  complex  demodulates  per  sequence.  The  complex 
demodulates  are  averaged  over  6  consecutive  sequences  in  1988  resulting  in  one  averaged 
sequence  record  every  184  seconds.  The  demodulates  are  averaged  over  4  sequences  in 
1989.  The  signal-to-noise  ratio  increases  by  101og106  =  8  dB  following  the  six-period 
average  in  1988  and  by  101og104  =  6  dB  following  the  four-period  average  in  1989. 
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D.2.2  Replica  Correlation 


One  period  of  a  (i(i.5  Hz  sampling  of  (0.1)  phase  angles  (expressed  as  complex  demodu¬ 
lates)  is  contained  in  the  following  sequence  of  length  J, 

*4>{j}  =  cos[|(rfch]  ~  ^)1  +  *sin[|(rfcb]  “  ^)]i  J  =  0,  1,...,  J  —  1,  (0.9) 

where, 

OO 

*b']  =  Z>b'-4/],  (0.10) 

i- 0 

and  u[j]  is  the  unit  step  sequence.  The  received  complex  demodulates  (D.6)  are  circularly 
convolved  with  a  sequence-reversed  (matched  filter)  replica  of  the  transmitted  demodu¬ 
lates  (0.9).  The  pulse-compressed  output  is  given  by, 

J- 1 

d[m]  =  m  =  0, 1, 1 .  (0.11) 

:=o 

where  the  periodic  sequence,  s^[j],  is  given  by, 

OO 

h[j]=  Y1  sAj  +  rJ}-  (0.12) 

r=— 00 

[Oppcnheim  and  Schafer ,  1989],  and  sj  is  the  complex  conjugate  of  $$■  The  replica 
correlation  (0.11)  compresses  30.7  s  of  energy  along  each  multipath  into  a  single  pulse 
of  one  digit  length  (~  0.06  s)  [Spiesberger  et  al.,  1989a],  whose  amplitude  and  location 
depend  on  the  gain,  bn,  and  delay,  rn,  in  (D.3).  The  signal-to-noise  ratio  increases  by 
1 0  log10  511  =  27  dB  following  pulse  compression,  and  a  sliding  average  of  4  complex 
samples  (separated  by  0.015  s)  increases  the  signal-to-noise  ratio  by  another  6  dB  at  the 
cost  of  some  resolution  [Spiesberger  et  al.,  1989a], 

D.2.3  Phase  Change  Between  Sequence  Records 

The  change  in  phase  of  the  mth  pulse-compressed  demodulate  from  record  r  -  1  to  r  is 
given  by  the  phase  of, 

d[m,  r  -  l]ri*[7/i,  r]. 

The  average  change  in  phase  from  record  r  -  1  to  r  is  given  by  the  ph 

M 

Ac'6'*’ lrl  =  d[m,r  -  l]<f*[m.r],  (D.13) 

rn =0 
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where  there  are  M  demodulates  containing  a  signal  [Spiesberger  et  ai,  1989a].  The  phase 
in  ( 1).  1 3 )  is  referred  to  as  the  cross-correlated  phase,  because  it  is  the  phase  of  the  complex 
cross  correlation  of  adjacent  pulse- com  pressed  records  at  zero-lag.  The  cross-correlated 
phase  can  be  converted  to  the  average  travel-time  difference  between  records  using  (2.1). 
The  precision  of  the  resulting  estimation  of  travel-time  difference  is  about  135  fis  [Spies- 
berger  et  ai,  1989a]. 
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